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Abstract 

We show in this paper that it is possible to attain very high, including observable, values 
for the level of non-gaussianity f^L associated with the bispectrum of the primordial 
curvature perturbation f, in a subclass of small-field slow-roll models of inflation with 
canonical kinetic terms. Such a result is obtained by taking care of loop corrections both 
in the spectrum P^ and the bispectrum B^. Sizeable values for /jvx arise even if Q is 
generated during inflation. Five issues are considered when constraining the available 
parameter space: 1. we must ensure that we are in a perturbative regime so that the £ 
series expansion, and its truncation, are valid. 2. we must apply the correct condition 
for the (possible) loop dominance in B^ and/or P^. 3. we must satisfy the spectrum 
normalisation condition. 4. we must satisfy the spectral tilt constraint. 5. we must have 
enough inflation to solve the horizon problem. 



1 Introduction 



Since COBE pQ discovered and mapped the anisotropics in the temperature of the cosmic 
microwave background radiation [2], many balloon and satellite experiments have refined the 
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measurements of such anisotropies, reaching up to now an amazing combined precision. The 
COBE sequel has continued with the WMAP satellite [3] which has been able to measure the 
temperature angular power spectrum up to the third peak with unprecedent precision [3], and 
increase the level of sensitivity to primordial non-gaussianity in the bispectrum by two orders of 
magnitude compared to COBE [HIE]. The next-to-WMAP satellite, PLANCK [7], whose launch 
is programmed for October 2008, is expected to precisely measure the temperature angular 
power spectrum up to the eighth peak [8], and improve the level of sensitivity to primordial 
non-gaussianity in the bispectrum by one order of magnitude compared to WMAP [9]. 

Because of the progressive improvement in the accuracy of the satellite measurements de- 
scribed above, it is pertinent to study cosmological inflationary models that generate significant 
(and observable) levels of non-gaussianity. An interesting way to address the problem involves 
the 5N formalism [10], [TT|, [12], which can be employed to give the levels of non-gaussianity f^L 
[T3] and t nl [HI [15] in the bispectrum and trispectrum of the primordial curvature per- 
turbation C, respectively. Such non-gaussianity levels are given, for slow-roll inflationary models, 
in terms of the local evolution of the universe under consideration, as well as of the n-point 
correlators, evaluated a few Hubble times after horizon exit, of the perturbations 5<pi in the 
scalar fields that determine the dynamics of such a universe during inflation. 

In the SN formalism for slow-roll inflationary models, the primordial curvature perturbation 
C(x, t) is written as a Taylor series in the scalar field perturbations 5<j)i(x, t+), evaluated a few 
Hubble times after horizon exit, 

C(t,x) = J2^ t ) s Mx,t*)-Y, N M(s<f>i(x,t*)) + 

i i 

ij ij 

E ^ i jfe(*)50i(x,t*)50j(x,t*)50 A .(x,t*) - — J2 N ijk (t){5(j) i (y.,U)5(j) j (yi,U)5(j) k (yi,U)) + 

ijk ' ijk 

+ ••• , (1) 

where the brackets mean spatial averages, N is the amount of inflation (or number of e-folds) 
from a bit later (in Hubble times) than the time when the cosmologically relevant scales exit the 
horizon and until the time at which one wishes to calculate £, and iVj = iV^ = J^.g^ , and so 
on. It is in this way that the correlation functions of ( (for instance, (CkiCk 2 Ck 3 )) can be obtained 
in terms of series, as often happens in Quantum Field Theory where the probability amplitude is 
a series whose possible truncation at any desired order is determined by the coupling constants 
of the theory. A highly relevant question is that of whether the series for 5N converges in 
cosmological perturbation theory and whether it is possible in addition to find some quantities 
that determine the possible truncation of the series, which in this sense would be analogous to 
the coupling constants in Quantum Field Theory. In general such quantities will depend on 
the specific inflationary model; the series then cannot be simply truncated at some order until 
one is sure that it does indeed converge, and besides, one has to be careful not to forget any 
term that may be leading in the series even if it is of higher order in the coupling constant. 
This issue has not been investigated in the present literature, and generally the series has been 
truncated to second- or third-order neglecting in addition terms that could be the leading ones 

[m H31 nn usi hhi lizi nn na 1201 isn [221 [23]. 
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The most studied and popular inflationary models nowadays are those of the slow-roll vari- 
ety with canonical kinetic terms [2U [25], [26] , because of their simplicity and because they easily 
satisfy the spectral index requirements for the generation of large-scale structures. One of the 
usual predictions from inflation and the theory of cosmological perturbations is that the levels 
of non-gaussianity in the primordial perturbations are expected to be unobservably small when 
considering this class of models [HI HH HH1 EQl [221 12ZI I2H [29l [30l [310. This fact leads us to 
analyze the cosmological perturbations in the framework of first-order cosmological perturbation 
theory. Non-gaussian characteristics are then suppressed since the non-linearities in the inflaton 
potential and in the metric perturbations are not taken into account. The non-gaussian char- 
acteristics are actually present and they are made explicit if second-order |34j or higher-order 
corrections are considered. 

The whole literature that encompasses the slow-roll inflationary models with canonical ki- 
netic terms reports that the non-gaussianity level f^L is expected to be very small, being of the 
order of the slow-roll parameters q and r/j, (e^, \r\i\ <C 1) [321 d21 EH1 E3 EE]- These works have 
not taken into account either the convergence of the series for ( nor the possibility that loop 
corrections dominate over the tree level ones in the n-point correlators. Our main result in this 
paper is the recognition of the possible convergence of the ( series, and the existence of some 
"coupling constants" that determine the possible truncation of the ( series at any desired order. 
When this situation is encountered in a subclass of small-field slow-roll inflationary models with 
canonical kinetic terms, the one-loop corrections may dominate the series when calculating ei- 
ther the spectrum P^, or the bispectrum B^. This in turn may generate sizeable and observable 
levels of non-gaussianity in total contrast with the general claims found in the present literature. 

The layout of the paper is the following: in Section [2] we consider the quantities that describe 
the statistical properties encoded in any probability distribution function; theoretical explana- 
tions as well as observational constraints for ( are given. Section [3] is devoted to the issue of the 
£ series convergence and loop corrections in the framework of the 5N formalism, as well as to the 
presentation of the current knowledge about primordial non-gaussianity in slow-roll inflation- 
ary models. A particular subclass of small-field slow-roll inflationary models is the subject of 
Section H] as it is this subclass of models that generate significant levels of non-gaussianity. The 
available parameter space for this subclass of models is constrained in Section by taking into 
account some observational requirements such as the COBE normalisation, the scalar spectral 
tilt, and the minimal amount of inflation. Another requirement of methodological nature, the 
possible tree-level or one- loop dominance in P^ and/or B^, is considered in this section. The 
level of non-gaussianity /tvl in the bispectrum B^ is calculated in Section [6] for models where 
( is generated during inflation; a comparison with the current literature is made. Section [7] is 
devoted to central issues in the consistency of the approach followed such as satisfying necessary 
conditions for the convergence of the ( series and working in a perturbative regime. Finally 

1 One possible exception is the two-field slow-roll model analyzed in Ref. [T7] (see also Refs. [35J[33]) where 
observable, of order one, values for /jvl are generated for a reduced window parameter associated with the initial 
field values when taking into account only the tree-level terms in both and Bq. However, such a result seems to 
be incompatible with the general expectation, proved in Ref. [18] . of f^L being of order the slow-roll parameters, 
and in consequence unobservable, for two-field slow-roll models with separable potential when considering only 
the tree-level terms both in P{ and B^. The origin of the discrepancy could be understood by conjecturing that 
the trajectory in field space, for the models in Refs. [171 [321 [33] , seems to be sharply curved, being quite near a 
saddle point; such a condition is required, according to Ref. [18], to generate fNL ~ 
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in Section [8] we conclude. The professional reader who is already familiarized with the present 
ideas on the cosmological non-gaussianity may skip Sections |2] and [31 leaping directly to the new 
material starting from Section HI As regards the level of non-gaussianity t^l in the trispectrum 
Tq, it will be studied following the sequence of ideas presented above in a companion paper [35J. 

2 Statistical descriptors for a probability distribution func- 
tion 

The primordial curvature perturbation C, as well as the contrast in the temperature of the cos- 
mic microwave background radiation ST/T and the gravitational potential $ 9 , are examples 
of cosmological functions of space and time being described by probability distribution func- 
tions. In particular, the probability distribution function /(C) for ( has well defined statistical 
descriptors which depend directly upon the particular inflationary model and that are suitable 
for comparison with present observational data. Such a comparison allows us either to reject 
or to keep particular inflationary models as those which better represent nature's behaviour. In 
this section we present a cosmologically motivated description of the statistical descriptors for 
probability distribution functions, focusing mainly on /(C)- 

2.1 Theoretical 

A probability distribution function /(C) for any function of space and time C(x, t) may be under- 
stood as the univocal correspondence between the possible values that C m ay take throughout 
the space and the normalised frecuency of appearences of such values for a given time. Any 
continous function of C might represent a probability distribution function as long as /(C) > 
and /^o/(C)^C = 1- However, for a particular probability distribution function, how many 
independent parameters do we need to completely characterize it in a unique way? And de- 
spite the possible infinite number of parameters required to do this, what is the information 
encoded in those parameters? The answers to these questions rely on the moments m^(n) of 
the distribution. 

For a given probability distribution function /(C), there are an infinite number of moments 
that work as statistical descriptors of C(x, t): 




the skewness 



the variance 



the kurtosis 








(3) 





and so on. 
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What can we say about C( x ;^) from the knowledge of the moments of the distribution? If, 
for instance, all the odd moments with n > 3 (skewness, ... etc) are zero, we can say that the 
probability distribution function f(() is even around the mean value. If in addition all the even 
moments with n > 4 (kurtosis, ... etc) are expressed only as products of the variance, we can say 
that the distribution function is gaussian. Indeed, as is well known, the only quantities required 
to reproduce a gaussian function are the mean value and the variance: 

f . (f\ = 1 p -(C-m c (l)) 2 /2m c (2) / fi x 

J gaussian VS J — / ° ■ V u / 

w27rm^(2) 

Departures from the exact gaussianity come either from no n- vanishing odd moments with n > 3, 
in which case the probability distribution function is non-symmetric around the mean value, or 
from higher even moments different to products of the variance, in which case the probability 
distribution function continues to be symmetric around the mean value although it is non- 
gaussian, or from both of them. A non-gaussian probability distribution function requires then 
more moments, other than the mean value and the variance, to be completely reconstructed. 
Such a reconstruction process is described for instance in Ref. [36] . 

Working in momentum space is especially useful in cosmology because the modes associated 
with the quantum fluctuations of scalar fields during inflation become classical once they leave 
the horizon [2l"l [371 [38] . The same applies for the primordial curvature perturbation ( which, 
in addition, is a conserved quantity while staying outside the horizon if the adiabatic condition 
is satisfied [12J. As regards the moments of the probability distribution function, they have 
a direct connection with the correlation functions for the Fourier modes Ck = / d 3 k((x)e _lk ' x 
defined in flat space. As the n-point correlators of Ck are generically defined in terms of spectral 
functions of the wavevectors involved^]: 

two — point correlation — > spectrum : 

(C kl Ck 2 ) = (27r) 3 5 3 (k 1 + k 2 )P c (A;), (7) 

three — point correlation — > bispectrum : 

(C kl Ck 2 Ck 3 > = (27rf8 3 (k 1 + k 2 + k 3 )B c (k 1 ,k 2 ,k 3 ), (8) 

four — point correlation — > trispectrum : 

(CkAsCkaCO = (2 7 r) 3 5 3 (k 1 + k 2 + k 3 + k 4 )T c (k 1 , k 2 , k 3 , k 4 ) , (9) 



and so on 



the moments of the distribution are then written as momentum integrals of the spectral functions 
for the modes Ck : 

(2 )3 ^(^ ' 



2 The homogeneity and isotropy requirements at large scales imply that the spectrum and bispectrum 
are functions of the wavenumbers only. For the trispectrum and the other higher order spectral functions, 
the momentum dependence also involves the direction of the wavevectors. 
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the skewness : m^(3) = / — TTTXe — B^(k\, k 2 , k 3 ) , (11) 



the kurtosis : m c (4) = / ~ T c( k i> k 2» k 3, k 4 ) , (12) 



(2tt) 9 



and so on. 

Non-gaussianity in ( is, therefore, associated with non-vanishing higher order spectral functions, 
starting from the bispectrum B^. 

Theoretical cosmologists work with (. However, astronomers work with observable quantities 
such as the contrast in the temperature of the cosmic microwave background radiation ST/T. 
The connection between the theoretical cosmologist quantity C and the astronomer quantity 
ST/T is given by the Sachs- Wolfe effect [39] which, at first-order and for superhorizon scales, 
looks as follows: 

= 4&- (13) 

Thus, although it is essential to study the Sachs- Wolfe relation at higher orders, which is far 
more complicated than Eq. ffl3|) . theoretical cosmologists may study the statistical properties 
of the observed 5T/T through the spectral functions associated with the curvature perturbation 

e 

mean value of 5T/T = — ► mean value of ( = , (14) 

variance : ms T / T (2) — > spectrum : P^(k) , (15) 

skewness : msT/T^) bispectrum: B^(k\, k 2 , ks) , (16) 

kurtosis: ms T / T (A) — > trispectrum : T^(ki, k 2 , k 3 , k 4 ) , (17) 




and so on. 

To end this section, we will parametrize the spectral functions of C in terms of quantities 
which are the ones for which observational bounds are given. Because of the direct connection 
between these quantities and the moments of the probability distribution function /(C), we may 
also call these quantities as the statistical descriptors for /(C)- The spectrum P^ is parametrized 

1 /2 

in terms of an amplitude and a spectral index which measures the deviation from an 
exactly scale-invariant spectrum [21]: 

™ - i^T" ■ (i8) 

where a is the global expansion parameter and H = a/ a is the Hubble parameter, with the 
dot meaning a derivative with respect to cosmic time. The bispectrum and trispectrum 



6 



are parametrized in terms of products of the spectrum P<^, and the quantities /jvl and t nl 
respectively [1HI27]: 



S c (A;i, k 2 , h) = -f NL [P c (fci)P f (& 2 ) + cyclic permutations] , (19) 

T^(ki, k 2 , k 3 , k 4 ) = ^tjvl [P f (/ci)P c (/c 2 )P^(|ki + k 4 |) + cyclic permutations] . (20) 

Higher order spectral functions would be parametrized in an analogous way. Given the present 
observational state-of-the-art, n^, /jvl, and tjvl are the statistical descriptors that discrimi- 

1/2 

nate among models for the origin of the large-scale structure once has been fixed to the 
observed value. Since non-vanishing higher order spectral functions such as and imply 
non-gaussianity in the primordial curvature perturbation (, the statistical descriptors /jvl and 
tnl are usually called the levels of non-gaussianity. 

2.2 Observational 

COBE provided us with a reliable value for the spectral amplitude V^ 2 |41j : V\j 2 = (4.957 ± 
0.094) x 10 -5 which is usually called the COBE normalisation. As regards the spectral index, 
the latest data release and analysis from the WMAP satellite shows that = 0.960 ± 0.014 [6] 
which rejects exact scale invariance at more than 2a. Such a result has been extensively used 
to constrain inflation model building [42], and although several classes of inflationary models 
have been ruled out through the spectral index, lots of models are still allowed; that is why it 
is so important an appropiate knowledge of the statistical descriptors /jvl and tjvl- Present 
observations show that the primordial curvature perturbation £ is almost, but not completely, 
gaussian. The level of non-gaussianity /jvl in the bispectrum B^, after five years of data from 
NASA's WMAP satellite, is in the range -9 < f NL < 111 at 2a [6J. There is at present 
no observational bound on the level of non-gaussianity tjvl in the trispectrum although 
it was predicted that COBE should either detect it or impose the lower bound \t nl \ < 10 8 
[TU 03] . It is expected that future WMAP data releases will either detect non-gaussianity or 
reduce the bounds on /jvl and tjvl at the 2a level to |/jvl| < 40 [9] and \tnl\ < 2 x 10 4 [H] 
respectively. The ESA's PLANCK satellite E], which will be launched in October 2008, 
promises to reduce the bounds to |/jvl| < 10 [9] and \tnl\ < 560 j II] at the 2a level if non- 
gaussianity is not detected. In addition, by studying the 21-cm emission spectral line in the 
cosmic neutral Hydrogen prior to the era of reionization, it is also possible to know about the 
levels of non-gaussianity /jvl and tjvl; the 21-cm background anisotropies capture information 
about the primordial non-gaussianity better than any high resolution map of cosmic microwave 
background radiation: an experiment like this could reduce the bounds on the non-gaussianity 
levels to | /jvl | < 0.2 [151 HE] and |tjvl| < 20 [IS] at the 2a confidence. Finally, it is worth stating 
that there have been recent claims about the detection of non-gaussianity in the bispectrum Bq 
of £ from the WMAP 3- year data [17J HH]. Such claims, which report a rejection of /jvl = at 
more that 2a (26.9 < /jvl < 146.7), are based on the estimation of the bispectrum while using 

3 There is actually a sign difference between the /jvl defined here and that defined in Ref. [57]. The origin of 
the sign difference lies in the way the observed /jvl is defined [S] , through the Bardeen's curvature perturbation: 
$ s = <f>£ + /atl(^l) 2 witn * B = ( 3 / 5 )C an( i the way f NL is defined in Ref. [57], through the gauge invariant 
Newtonian potential: <J> W = + /jvl($l f with $W = ~( 3 / 5 )C US- 
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some specific foreground masks. The WMAP 5- year analysis [6] shows a similar behaviour when 
using those masks, but reduces the significance of the results when other more conservative 
masks are included allowing again the possibility of exact gaussianity. 



3 The SN formalism 

The SN formalism [TUl [TT1 fl~2] provides a powerful method for calculating £ and all its statistical 
descriptors at any desired order in cosmological perturbation theory from knowing, in the case 
of slow-roll inflationary models, only the evolution of a family of unperturbed universes and 
the correlation functions a bit later than horizon exit of the perturbations in the field scalars 
present during inflation [13] . This section is devoted to a brief review of the formalism and to a 
short discussion of some relevant issues which either have not been properly taken into account 
or have not been discussed at all in the current literature. 



3.1 The non-linear primordial curvature perturbation ( 



In the cosmological standard model [19] the observable Universe is homogeneous and isotropic, 
being described by the unperturbed Friedmann- Robertson- Walker metric whose line element, 
for a spatially flat universe, looks as follows: 



ds 2 = -dt 2 + a 2 (t)S ij dx i dx j 



(21) 



where a(t) is the global expansion parameter, t is the cosmic time, and x represents the position 
in cartesian spatial coordinates. The homogeneity and isotropy conditions describe very well the 
Universe at large scales, but departures from the unperturbed background are observationally 
evident at smaller scales. 

One way to parametrize the departures from the homogeneous and isotropic background is 
to include perturbations in the metric, for which we have to define a slicing and a threading. 
The slicing will be defined so that the energy density in fixed-t slices of spacetime is uniform. 
The threading will correspond to comoving fixed-x world lines. Thus, the perturbed spatial 
metric may be defined as 

gy = a*{t)tXto*hv(t,x) , (22) 
where 7y(i, x), which gives origin to the tensor perturbations, has unit determinant. This means 
that we may define a local scale factor 



a(t, 



a(t)e 



C(*,x) 



(23) 



which is proportial to any volume of the Universe bounded by fixed spatial coordinates. 

The interesing feature of Eq. (1231) is that it is possible to write the £(t, x) quantity in terms 
of the perturbation in the amount of expansion from an initial time t in where the slice is flat 
(i.e. the spatial metric in such a slice is the same as that in Eq. fl22l) but without the £ factor) 
to a final time t where the slice is of uniform energy density: 



((t,xL) = 6N = N{t,x) - N (t) 



In 



a(t, x) 
. a(t in ) 



In 



' a{t) 
a(t in ) 



(24) 
(25) 
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This is the 5N formalism (101 HU [12] where the ((t, x) quantity, being a non-linear extension of 
the primordial curvature pertubation, reproduces the usually defined gauge-invariant curvature 
perturbation (i at first order [50] : 

Ci = -V - ^— , (26) 
P 

where p is the energy density and ip is the scalar perturbation in the spatial metric at super- 
horizon scales and at first order: 

g ij = a 2 {t){l-2^)8 ij . (27) 

(, as non-linearly defined in Eqs. (j22p and ([Ml), is a conserved quantity on superhorizon scales 
as long as the adiabatic condition (the pressure as a function of only the energy density) is 
satisfied [12]. 



3.2 £ series convergence and loop corrections 

In order to calculate C(^, x ) from Eq. f[24l) . we need information about the physical content of 
the Universe at times t and t- m . By choosing the initial time t- m a few Hubble times after the 
cosmo logically relevant scales leave the horizon during inflation t m = t*, and the final time t 
corresponding to a slice of uniform energy density, we recognize that N, for slow-roll inflationary 
models, is completely parametrized by the values a few Hubble times after horizon exit of the 
scalar fields fa present during inflation and the energy density at the time at which one wishes 
to calculate (: 

C(t, x) = N{p{t), fa%, x), fa(t+, x), ...) - N{p{t), fa(t*), <h(U), •••) • (28) 

The previous expression can be Taylor-expanded around the unperturbed background values for 
the scalar fields fa and suitably redefined so that x)) = 0. Thus, 

C(t,x) = J2 N M s Mt*,x)-T, N M( s Mt*,*)) + 

i i 

+^£%(*)W*> x )^(**> x ) -|e%(*)<^(**> x )%(^ x )> + 

ij ij 

+^ E Nij k (t)5fa(t±, x)5(f)j(U, x)84> k (t*, x) - ^ ^ N ijk(t)(Sfa(U, x) x)50 fc (i*, x 

ijk ' ijk 

+ ... , 

where the 5</>j(t*,x) are the scalar field perturbations in the flat slice a few Hubble times after 
horizon exit, whose spectrum amplitude is given by 



< = § . (30) 



and the notation for the N derivatives is iVj = J^, iVy = J^.q^. , and so on. 
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The expression in Eq. fl29l) has been used to calculate the statistical descriptors of ( at any 
desired order in cosmological perturbation theory by consistently truncating the series [13]. For 
instance, by truncating the series at first order, the amplitude of the spectrum of ( defined 
in Eqs. ([7J) and <^ is given by p] 



n = (§) 2 E^, (31) 



which in turn gives the well known formula for the spectral index |llj : 

9 Em ViNjNa . , 

where a subindex i in V means a derivative with respect to the <f>i field, and being e one of the 
slow-roll parameters defined by e = —H/H 2 , mp = (8ttG)~ 2 the reduced Planck mass, and V 
the scalar inflationary potential. Analogously, the level of non-gaussianity /jvx i n the bispectrum 
of ( defined in Eqs. (EJ) and (fl9l) is obtained by truncating the series at second order and 
assuming that the scalar field perturbations 5(j)i are perfectly gaussian [13] : 



6 Eij N,NjNjj j: m N,,N, k N kl 

E,^f +n E.A5T M * £) ' (33) 

In the last expression the ln(/cL) factor is of order one, L being the infrared cutoff when calcu- 
lating the stochastic properties in a minimal box [521 [53] . 

The truncated series methodology has proved to be powerful and reliable at reproducing 
successfully the level of non-gaussianity /nl in single-field slow- roll models [29J and in the 
curvaton scenario [13] . Nevertheless, for more general models, how reliable is it to truncate the 
series at some order? In the first place, from Eq. (1291 it is impossible to know whether the 
series converges until the N derivatives are explicitly calculated and the convergence radius is 
obtained; obviously if the series is not convergent at all, the expansion in Eq. (129]) is meaningless. 
Without any proof of the contrary, the current assumption in the literature [HI H3J HH [15j [16j 
ITT] [181 [T9| 1201 I2TI 122] [23] has been that the ( series is convergent. In addition, supposing that 
the convergence radius is finally known, the truncation at any desired order would again be 
meaningless if some leading terms in the series get excluded. Such a situation might easily 
happen if each x-dependent term in the ( series is considered smaller than the previous one, 
which indeed is the standard assumption [IH [TBI H31 H51 USl HZl HSl HSl HOI I2H 1221 1^] , but which 
is not a universal fact. 

When studying the series through a diagrammatic approach [51] , in an analogous way to that 
for Quantum Field Theory via Feynman diagrams, the first-order terms in the spectral functions 
are called the tree-level terms. Examples of these tree-level terms are those in Eqs. (j3T|) and 
([32]) . and the first one in Eq. (1331) . Higher-order corrections, such as that which contributes with 
the second term in Eq. ([33]) . are called the loop terms because they involve internal momentum 
integrations. The statistical descriptors of ( has been so far studied by naively neglecting the loop 
corrections against the tree- level terms [TIllSllIllliailSllIZllIBllISllSnilSIl^l^; nevertheless, 
as might happen in Quantum Field Theory, eventually some loop corrections could be bigger 
than the tree-level terms, so it is essential to properly study the possible n-loop dominance in 
the spectral functions. 
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3.3 Non-gaussianity in slow-roll inflation 

The most frequent class of inflationary models found in the literature are those which satisfy 
the so called slow-roll conditions, as these very simple models easily meet the spectral index ob- 
servational requirements discussed in Subsection 12.21 for the generation of large-scale structures. 
The slow-roll conditions for single-field inflationary models with canonical kinetic terms read 

4> 2 « V(<j>), (34) 
101 « W|, (35) 

where is the inflaton field and V(<p) is the scalar field potential. On defining the slow-roll 
parameters e and rj^ as [21] 

' = -Jj. (36) 

* s f --w (37) 

the slow-roll conditions in Eqs. (1341) and fl35|) translate into strong constraints for the slow-roll 
parameters: e, \t)^,\ <C 1, which actually become e, \r}^\ < 10~ 2 in view of Eq. fl32l for single-field 
inflation: 

n ( -l = 2^-6e, (38) 

and the observational requirements presented in Subsection 12 .21 

Multifield slow-roll models may also be characterized by a set of slow-roll parameters which 
generalize those in Eqs. ( 136]) and ( 1371) [25J: 

m = ™\ V f- (40) 

By writing the slow-roll parameters in terms of derivatives of the scalar potential, as in the last 
two expressions, we realize that the slow- roll conditions require very flat potentials to be met. 

The level of non-gaussianity J'nl in slow-roll inflationary models with canonical kinetic terms 
has been studied both for single- field [27] and for multiple- field inflation [181 UHl [20J , assuming ( 
series convergence and considering only the tree- level terms both in and in B^. What these 
works find is a strong dependence on the slow-roll parameters 6j and rji] for instance, Ref. [27] 
gives us for single-field models: 

6 

-f NL = e(l + f) + 2e- V(P , (41) 

where / is a function of the shape of the wavevectors triangle within the range < / < 5/6. 
Refs. [221 ES] show that in such a case the inclusion of loop corrections is unnecessary because 
the latter are so small compared to the tree-level terms. Thus, f^L in single-field models with 
canonical kinetic terms is slow-roll suppressed and, therefore, unobservably small. As regards the 
multifield models, J'nl was shown, first in the case of two-field inflation with separable potential 
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[T8] and later for multiple-field inflation with separable [19] and non-separable (20] potentials, 
to be a rather complex function of the slow-roll parameters and the scalar potential that in 
most of the cases ends up being slow-roll suppressed. Only for models with a sharply curved 
trajectory in field space might the /nl be at most of order one, the only possible examples to 
date being the models of Refs. [T71 [32], [33] (see anyway footnoted]). Again, such predictions are 
based on the assumptions that the ( series is convergent and that the tree-level terms are the 
leading ones, so they might be badly violated if loop corrections are considered. 

Following a treatment parallel to that in Ref. [IB] , the level of non-gaussianity t^l is cal- 
culated in Ref. [22] for multifield slow-roll inflationary models with canonical kinetic terms, 
separable potential, and assuming convergence of the ( series and tree-level dominance. From 
reaching similar conclusions to those found for the f^L case, the t nl is slow-roll suppressed in 
most of the cases although it might be of order one if the trajectory in field space is sharply 
curved. Nevertheless, as will be shown in a companion paper [35] , there may be a big enhance- 
ment in tnl if loop corrections are taken into account. 

Finally, it is worth mentioning that there are other classes of models where the levels of 
non-gaussianity /tvl and r NL are big enough to be observable. Some of these models correspond 
to general langrangians with non-canonical kinetic terms (k- inflation [56], DBI inflation |57j . 
ghost inflation [58j, etc.), where the sizeable levels of non-gaussianity have mostly a quantum 
origin, i.e. their origin relies on the quantum correlators of the field perturbations a few Hubble 
times after horizon exit. Non-gaussianity in has been studied in these models for the single- 
field case [59], ED] and also for the multifield case [HH E2] [63], E3]- A recent paper discusses 
the non-gaussianity in Tq for these general models for single-field inflation [65]. In contrast, 
there are some other models where the large non-gaussianities have their origin in the field 
dynamics at the end of inflation [66l [67] ; nice examples of this proposal are studied for instance 
in Refs. [6H1 E3 E2 EE]. However, since the inflationary models of the slow-roll variety with 
canonical kinetic terms are the simplest, the most popular, and the best studied so far although, 
in principle, the non-gaussianity statistical descriptors are too small to ever be observable, it 
is very interesting to consider the possibility of having an example of such models which does 
generate sizeable and observable values for f^L- This appealing possibility will be the subject 
of the following sections. 

4 A subclass of small-field slow-roll inflationary models 

According to the classification of inflationary models proposed in Ref. [72], the small- field models 
are those of the form that would be expected as a result of spontaneous symmetry breaking, 
with a field initially near an unstable equilibrium point (usually taken to be at the origin) and 
rolling toward a stable minimum ((f)) ^ 0. Thus, inflation occurs when the field is small relative 
to its expectation value (f> <C ((f)). Some interesting examples are the original models of new 
inflation [73], [73], modular inflation from string theory [75], natural inflation [76], and hilltop 
inflation [77] . As a result, the inflationary potential for small-field models may be taken as 




(42) 
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Figure 1: Our small- field slow- roll potential of Eq. (j43p with rj^, r\ a < 0. The inflaton starts near 
the maximum and moves away from the origin following the o = trajectory depicted with the 
solid black line. (This figure has been taken from Ref. [T7]). 



where the subscript i here denotes the relevant quantities of the ith. field, p is the same for all 
fields, and Aj and fii are the parameters describing the height and tilt of the potential of the ith 
field. 

While Ref. [TS] studies the spectrum of ( for general values of the parameter p and an 
arbitrary number of fields, assuming ( series convergence and tree-level dominance, we will 
specialize to the p = 2 case for two fields and a: 

v = v "( 1 + b*& P + l"^)- (43) 

where we have traded the expressions 

A x +A 2 for V , (44) 



and 



, for -Vq^V, (45) 
/if 2m z p 



H for -V -\. (46) 
H2 zrrip 

On doing this, and assuming that the first term in Eq. (j43p dominates, rj^ < and r/ CT < 
become the usual rj slow-roll parameters associated with the fields <j) and <J. 
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We have chosen for simplicity the a = trajectory (see Fig. [T]) since in that case the potential 
in Eq. fj4*3|) reproduces for some number of e-folds the hybrid inflation scenario [79] where <ft is 
the inflaton and a is the waterfall field. Non-gaussianity in such a model has been studied in 
Refs. [131 E3 Ed EU EOj El]; in particular, Ref. [13] used a one-loop correction to conjecture 
that Jnl in this model would be sizeable only if ( was not generated during inflation, which 
turns out not to be a necessary requirement as we will show later. Ref. [TTJ, in contrast, works 
only at tree-level with the same potential as Eq. (|4*3|) but relaxing the a = condition, finding 
that values for f NL ~ Oil) are possible for a small set of initial conditions and assuming a 
saddle-point like form for the potential < and r\ a > 0). 



5 Constraints for having a reliable parameter space 

We will explore now the constraints that the model must satisfy before we calculate Jnl- Our 
guiding idea will be the consideration of the role that the tree-level terms and one-loop cor- 
rections to both and have in the determination of the available parameter space. Only 
after calculating f NL in Section [6] will we come back to the discussion of the consistency of the 
approach followed in the present section by studying the ( series convergence and the validity 
of the truncation at one loop level. 

5.1 Tree-level or one- loop dominance 

Since we are considering a slow-roll regime, the evolution of the background and a fields in 
such a case is given by the Klein-Gordon equation 

4> + 3#0 + Vj, = , (47) 

supplemented with the slow-roll condition in Eq. (|35l) . This leads to 

0(iV) = ^exp(-iV^), (48) 
a(N) = a*ex P (-N Va ), (49) 

so the potential above leads to the following derivatives of N with respect to 0* and a* for the 
a = trajectory^ 

Ntj, — — ^— , N a = , (50) 



N H = -— - , N 4Hr = 0, N aa = -g- exp[2iV( % - rj a )] , (51) 



2 2r? 2 
N 4+t>=ZnH> N M<t = °> N au(j) = — g^exp[2JV(^-^)], N aaa = , (52) 



and so on. 



4 



When calculating the iV-derivatives, we have considered that the final time corresponds to a slice of uniform 
energy density. This means, in the slow-roll approximation, that V is homogeneous. 
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By means of the 5N formalism, we can make use of the above formulae to calculate the 
spectrum and the bispectrum of the curvature perturbation including the tree-level and the 
one-loop contributions when j^l > \r}^\ (see Appendix [A}. This is the interesting case since, as 
will be shown in Section El it generates sizeable values for f^L- Following the results in Appendix 
lAl we will write down just the leading terms to the tree- level and one- loop contributions given 
in Eqs. flUD, AMD, flATLD, and ([A28]): 

£V ■ (53) 

(kL) , (54) 





1 




spl—loop 




4^ 


jjtree 


i 




jjl — loop 

B ( 


V* 









1 - 1^1)] 


(i) 



HA *4«*(M), (55) 



2W \Uikf 
exp[6iV(K| - |^|)] (^) 6 ln(fcL)4vr 4 (^^j . (56) 

Because of the exponential factors in Eqs. (154")) and (156]) . it might be possible that the one- 
loop corrections dominate over P^ and/or B^. There are three posibilities in complete connection 
with the position of the <fr held when the cosmologically relevant scales are exiting the horizon: 

5.1.1 Both and P^ are dominated by the one-loop corrections 

Comparing Eqs. (1551) with (1541) and Eqs. (155)) with (155]) we require in this case that 

4ex P [4A% CT |-|^|)] » / 2 , (57) 

1 / XI -A- 



7? 



2tt 



^exp^M-M)] » 1 2 . (-»> 

1 / XI* 



2tt 



in which case only the first inequality is required. Employing the definition for the tensor to 
scalar ratio r |26j : 

Vt = i%K**) (59) 

1 12 

V T being the amplitude of the spectrum for primordial gravitational waves, we can write such 
an inequality as 

<<^|exp[4iV(|r /(J |-M)]. (60) 
\mpj 8 i}^ 

From now on we will name the parameter window described by Eq. (1601) as the low 0* region 
since the latter represents a region of allowed values for </>* limited by an upper bound. 
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5.1.2 dominated by the one-loop correction and P^ dominated by the tree-level 
term 

Comparing Eqs. ( |53j) with ( 1541) and Eqs. ( 1551) with ( 156]) we require in this case that 



gexp[4iV(|^|-|^|)] « - ^ x2 , (61) 

2^ 



%exp[6JV(|^|-N)] » ; * 2 . ( 62 ) 
77; 1 f -H* 



2tt 



which combines to give, employing the definition for the tensor to scalar ratio r introduced in 
Eq. (J59D, 

^4^P[4iV(k| - |^|)] « « ^4 ex P[ 6Ar d^l - 1^1)]- ( 63 ) 

From now on we will name the parameter window described by Eq. ( 1631) as the intermediate 0* 
region since the latter represents a region of allowed values for 0* limited by both an upper and 
a lower bound. 



5.1.3 Both B( and P^ are dominated by the tree-level terms 

Comparing Eqs. fl53l) with fl54l) and Eqs. (155]) with (156]) we require in this case that 

4exp[4iy(|r7 CT |-M)] « - * 2 , (64) 
"* A (§) 

^expfeiVd^l-^i)] « / 2 , (65) 



7? 



2tt 



in which case only the second inequality is required. Employing the definition for the tensor to 
scalar ratio r introduced in Eq. fl59l) . we can write such an inequality as 

(—) »^4 ex P[ 6Ar (W-M)l- (66) 
\mpj 8 r% 

From now on we will name the parameter window described by Eq. f l6"6"|) as the high <fr+ region, 
since the latter represents a region of allowed values for 0* limited by a lower bound. 



5.2 Spectrum normalisation condition 

The model must satisfy the COBE normalisation on the spectrum amplitude [UJ consider- 
ing that ( is assumed in this paper to be generated during inflation^. There exist two possibilities 
discussed right below. 

5 The scenario where £ is assumed not to be generated during inflation will be presented in a companion paper 
[35]. 
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5.2.1 £ generated during inflation and P^ dominated by the one-loop correction 

According to Eq. (l54l) . and the tensor to scalar ratio r definition in Eq. ( l59l) . we have in this 



case 



V\- lo ° v = ^ex P [AN(\ Va \-\ V ,\)](^)\n(kL) 



9 / \ 4 sr* 2 

% r ,. . . m m P \ frVt^' 



Lf exp[4iV(|^| - |^|)] ( -Jl ) ( -Li ) \n(kL) , (67) 



which reduces to 



•i 



2 ^2 



' > P c ^exp[4iV(|^| - 1^1)] \n(kL) , (68) 



where must be replaced by the observed value presented in Subsection 12.21 

5.2.2 ( generated during inflation and P^ dominated by the tree-level term 

According to Eq. (1531) . and the tensor to scalar ratio r definition in Eq. (1591) . we have in this 
case 

Tjtree \ ( 

' " ~ ~ 27T 

2 

(69) 



C " V2 



1 / mp\ rV( 



which reduces to 

\m P ) 7^8 

Notice that in such a situation, the value of the (f> field when the cosmologically relevant scales 
are exiting the horizon depends exclusively on the tensor to scalar ratio r, once rj^ has been 
fixed by the spectral tilt constraint as we will see later. 

5.3 Spectral tilt constraint 

The combined 5-year WMAP + Type I Supernovae + Baryon Acoustic Oscillations data [6] 
constrain the value for the spectral tilt as 

n ( - 1 = -0.040 ±0.014. (71) 

Here again we have two possibilities: V( is dominated either by the one-loop correction or by 
the tree-level term: 
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5.3.1 dominated by the one-loop correction 

In this case the usual spectral index formula at tree-level [TT] gets modified to account for the 
leading one-loop correction: 

n < - 1 = - 4£ - 2m -lit + ml ■ (72) 

By making use of the derivatives in Eqs. (1501) . (loTj) . and (1521) . we have 

n c - 1 = -4e + 4r/ ff + [hi^L)]" 1 , (73) 

which implies that the observed value for is never reproduced in view of In(kL) ~ 0(1). 
Moreover, when calculating the running spectral index dn^/dlnk from Eq. (173|) . we obtain 

^ C - - [ln(fcL)]" 2 , (74) 



(ilnfc 



which rules out the possibility that P^ is dominated by the one-loop correction since the 
calculated dn^/dlnk is far from the observationally allowed 2a range of values: —0.0728 < 
dnjdlnk < 0.0087 [60. 

5.3.2 P^ dominated by the tree-level term 

Now the usual spectral index formula [11] applies: 

n f - 1 = -2e - 2m P ^ ^ , (75) 

giving the following result once the derivatives in Eqs. (J50|) . (15T1) . and (152|) have been used: 

n c - 1 = -2e + 2?7 . (76) 

The efect of the e parameter may be discarded in the previous expression since, as often happens 
in small-field models [HI [77], e is negligible being much less than \r] a \: 



6 2 V 2 



1 ' ' 



2 \mp 



< M < \Va\ , (77) 



according to the prescription that the potential in Eq. (1431) is dominated by the constant term. 
Thus, by using the central value for — 1, we get 

^ = -0.020 . (78) 



6 We thank Eiichiro Komatsu for pointing out to us the dependence of uq and dn^/dlnk on ln(fcL). 
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5.4 Amount of inflation 



It is well known that the number of e-folds of expansion from the time the cosmological scales 
exit the horizon to the end of inflation is presumably around but less than 62 [24"l [82] . The 
slow-roll evolution of the field in Eq. (I48p tells us that such an amount of inflation is given by 

N= J-\ n (^A < 6 2. (79) 



Because of the characteristics of the inflationary potential in Eq. (l4~3]) . there is no definite 
mechanism for ending inflation in this model. It could not be by means of the violation of the 
e < 1 condition since this would imply extrapolating our results to a region where the potential 
in Eq. (1431) is no longer dominated by the constant term which, in addition, would spoil the large 
non-gaussianity generated and could send the model to an unknowable quantum gravity regime. 
Keeping in mind the results of Ref. [83] which say that the ultraviolet cutoff in cosmological 
perturbation theory could be a few orders of magnitude bigger than mp, we will therefore assume 
that inflation comes to an end when |v7</,|0 2 jimp ~ 10~ 2 . This allows us to be on the safe side 
(avoiding large modifications to the potential coming from ultraviolet cutoff-suppressed non- 
renormalisable terms, and keeping the potential dominated by the constant Vo term), leaving 
the implementation of a mechanism for ending inflation for a future wori0. Coming back to Eq. 
(179|) . we get then 

(80) 

which leads to 

^ > exp(-62|^|). (81) 



mp 



6 /nl 

In this section we will calculate the level of non-gaussianity represented in the parameter f^L by 
taking into account the constraints presented in Subsections 15.21 15.31 and 15.41 and the different 
0* regions discussed in Subsection 15.11 



6.1 The low (j)+ region 

This case is of no observational interest because dominated by the one-loop correction is 
already ruled out by the observed spectral index and its running as shown in Subsubsection l5~3.ll 
In addition, the generated non-gaussianity is so big that it causes violation of the observational 

7 We hope that the implementation of such a mechanism in our model will keep, or perhaps enhance, the 
generated non-gaussianity. Nevertheless the opposite behaviour might as well happen. For instance, Ref. [84] 
studies within a stochastic formalism a quadratic two-component slow-roll model without a dominant constant 
term in the potential. A momentary violation of the slow-roll conditions around the end of inflation shows 
enhancement of f^L to observable levels; however, such an enhancement vanishes once inflation ends completely. 
These results have been confirmed numerically within the SN formalism in Refs. [181 121) . 
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constraint Jnl > — 9: 



\f NL = c 3 = -[V C ln(£L)]~ 1/2 ~ -2 x 10 4 , (82) 

1 Li i 

according to the expressions in Eqs. (|T8|) . (|19l . (154|) . and (|56|) . 

We want to remark that, although it is of no observational interest, this case represents the 
first example of large non-gaussianity in the bispectrum of ( for a slow-roll model of inflation 
with canonical kinetic terms. It is funny to realize that the model in this case is additionally 
ruled out because the observational constraint on Jnl is violated by an excess and not by a 
shortfall as is currently thought [TBI [191 [201 [271 [28] . 

6.2 The intermediate 0* region 

The level of non-gaussianity, according to the expressions in Eqs. (TT8"j) . (fl^|) . (151I1) . and (151)1) . is 
in this case given by 

TfNL = — ^ = ^-exp[6iV(|^|-M)](-^) \n(kL) 

= ^exp[6JV(|i^| - |^|)] ^ 2 ^ln(A;L) 

= ^exp[6JV(|^| - \ V4> \)}V c ln(kL) , (83) 

r/arL ~ -2.457 x 10" 9 |r7 ff | 3 exp[300 ln(5.657 x 10~V 1/2 ) (|^| - 0.020)] , (84) 



where in the last line we have used expressions from Eqs. flTUj) . (1781 . and f[8"Uj) . 

Now, by implementing the spectral tilt constraint in Eq. (ITS]) in the spectrum normalisation 
constraint in Eq. (f70l) and the amount of inflation constraint in Eq. (jHTT) . we conclude that the 
tensor to scalar ratio r is bounded from below: r > 2.680 x 10~ 4 . 

In the plot r vs in figure [2] we show lines of constant f^L corresponding to the values 
f NL = —5,-10,-15. We also show the high and intermediate 0* regions in agreement with the 
constraint in Eq. f}63l : 

^4exp[4A% CT | - M)] « « ^4 ex P[ 6Ar d^l - > 

=>• 8.139 x 10 6 < l^l 3 exp[300 ln(5.657 x 10" 2 r" 1/2 ) (|r/ ff | - 0.020)] < 8.210 x 10 12 . (85) 

As is evident from the plot, the WMAP (and also PLANCK) observationally allowed 2a range of 
values for negative /nl, —9 < /nl, is completely inside the intermediate 0* region as required. 
More negative values for Jnl, up to Jnl — —20.647 are consistent within our framework for the 
intermediate 0* region, but they are ruled out from observation. Nevertheless, like for the low 0* 
region studied above, it is interesting to see a slow-roll inflationary model with canonical kinetic 
terms where the observational restriction on Jnl may be violated by an excess and not by a 
shortfall. So we conclude that if is dominated by the one-loop correction but Pq is dominated 
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Figure 2: Contours of f^i in the r vs \r] a \ plot. The intermediate (high) 0* region is below 
(above) the LIf line. The WMAP (and also PLANCK) observationally allowed 2a range of 
values for negative /jvl, —9 < f^L, is completely inside the intermediate 0* region. Notice that 
the LIf line matches almost exactly the f^L = —1.667 line. 



by the tree-level term, sizeable non-gaussianity is generated even ifQ is generated during inflation. 
We also conclude, from looking at the small values that the tensor to scalar ratio r takes in figure 
[2] compared with the present technological bound r > 10~ 3 [85J, that for non-gaussianity to be 
observable in this model, primordial gravitational waves must be undetectable. 

Notice that in oder to get positive values for Jnl, which is observationally more interesting 



in view of the results presented in Subsection I2.2[ r\ a should be positive according to Eq. fl83|) . 
However, being rj^ negative in order to reproduce the observed spectral tilt, the argument of 
the exponential in Eq. (53} would be negative, making the /nl obtained too small to be 
observationally interesting^. As regards the general case, in view of the previous reason being 
model dependent, we may only say that in order to get f^L positive when Bq is dominated by 
the one-loop corrections, should be positive (based on the definition of f^L m Eq. (fT9l) ) 
which means that the maximum between N a(T and should be positive in view of Eq. (1A12j) . 



Finally we want to point out that, by reducing our model to the single-field case, the con- 
sistency relation between f NL and presented in Ref. |86j: f^L ~ C( n c — -0 ^ s no ^ violated 
since in that case Bq is never dominated by the one-loop corrections for slow-roll inflation as 
demonstrated in Ref. [S2jE Thus, the level of non-gaussianity /nl for our model reduced to the 



8 We thank Eiichiro Komatsu for questioning us about this issue. 
9 We thank Filippo Vernizzi for questioning us about this issue. 
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single-field case is described by the high 0* region as shown below. 



6.3 The high 0* region 

This case is of no observational interest because, according to the expressions in Eqs. ( fl8l) . ( fl9l) . 
( 1531) . ( 1551) . and ( 1781) . the non-gaussianity generated is too small to be observable: 

TfNL = yr~^ = "Tfo = 0.020 , (86) 

in agreement with the consistency relation of Ref. [86] for our model reduced to the single-field 
case, and with the general expectations of Refs. [HI HH [201 123 [28] for slow-roll inflationary 
models with canonical kinetic terms where only the tree-level contributions are considered. 



7 Convergence of the ( series and perturbative regime 

In Sections [5] and [6] we have worked up to the one-loop diagrams in order to constrain the 
parameter space and find the level of non-gaussianity f^L- It is time then to address the issue of 
the ( series convergence and justify the existence of a perturbative regime so that the truncation 
of the series up to the one-loop order, for the model we have considered, is valid. A way to do 
that is by rederiving the ( series in terms of 50* and 5cr* by equating the unperturbed scalar 
potential to the perturbed one at the final time t; this of course is valid in view of the first 
slow-roll condition in Eq. (IMI) and the final slice being one of uniform energy density: 



i w A-exp[-2^] + ir/ ff ^exp[-2iVr ?(7 ]) 
2 uip 2 mp J 



V 



1 + 



+ 



] - exp[-2(iV + SN)^] + l ^ + f 7 *) 2 exp[-2(iV + 5N) V<7 ] 



m. 



From the previous expression it follows that 

rj^l exp[-2A^?7 9i ] + r] a al exp[-2iV7fc] 

= 7^(0* + 50*) 2 exp[-2(iV + SN)^] + T] a {a* + 5a*) 2 exp[-2(iV + 5N)r] a ] 
which is easier to handle in terms of variables x and y defined as 

50, 



X = 



rf a a 2 



1/2 



i7) 



(88) 



(89) 



(90) 



Thus, the exponentials factors contaning N (but not SN) are completely absorbed in y, and the 
expression in Eq. ( 1HH1) looks as follows: 



1 + 



8a* 



1 + x) 2 exp[-25iV^] + !| 2/ 2 exp[-25iVr ?(T ] . 



(91) 
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If we were able to solve for SN in Eq. fl9Tj) in terms of rj^, rj a , x, and y (after making = 0), 
we could Taylor-expand around x = and y = reproducing the x-dependent part of Eq. fl29|) . 
This would be really good because the Taylor expansion would look so clean, in the sense that 
all the concerning exponential factors contaning N which appear explicitely in Eq. ( |29i) would 
already be absorbed in y, that the issue of truncating at some specific order in 5<f)+ and 5<j± 
would be simply justified by requiring \x\ <C 1 and \y\ <C 1. Nevertheless, as is seen in Eq. (|9T|) . 
it is impossible to solve for 5N in terms of T]^, r] a , x, and y unless we make a Taylor expansion 
of the exponential functions aroud SN = 0: 







[1 + xf 



, ^ 2 

+ -%y 

Tor 



1 + 



+SN 



-2^(1 + x) 



2^ 2 



SN 2 



+ 



2 V l(l+x) 2 + 2 V ly 2 



+ 



(92) 



Notice that the Taylor expansion of the exponential functions is always convergent whatever the 
arguments of the exponentials are. Moreover, if the Taylor expansion derived from a function 
f(x) converges, it converges precisely to f(x) [87J. Thus, the expression in Eq. (l92il is actually 
the same as the expression in Eq. fl9TT) . 

Now, solving for SN in terms of rj^, rj a , x, and y, although possible in view of the expression in 
Eq. (192]) . is not an easy business. That is why we will truncate the series in Eq. (1921) up to second 
order" SN and solve the resultant quadratie equatiofl Notiee that , since ( = SN — (SN) 
and ( ~ 10~ 5 , we may truncate the series in Eq. ( 1921) up to whatever order we wish and still 
reproduce ( to high accuracy. Thus, the solution for the quadratic equation coming from the 
series in Eq. (|92|) after truncation at second order is: 

i "l 2 

,2 , 2 



5N 




1 , \2 . 1 V(f> 2 

l + x) +-—y 
2rj a 



1 + x 



2^ 



1 + x 



l + x) 2 + y 2 



2 

i V<t> 2 
-1 



§£*)' 



:i+x) 



X 



(93) 



If in addition we make Taylor expansions of the square root and the factor in the third line of 
the previous expression around x = and y = 0: 

1/2 



l + x) 



l V<t> 2 

2r]a 



(l+x) 



i-x 2 + ^M 
2 2t7 ct 



2 

+ 



y 2 + 



l 



l 



;i + x) 2 + y 2 







'(l + x) 2 + / 







l-2x + 3x 2 -y 2 + ... 



(94) 



(95) 



10 The truncation up to second order in 8N has been chosen in order to have complete consistency with the 
order of the variables x and y in Eq. (pJTj) . 
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introducing them into Eq. fl93|) . we end up with the following power series for 5N: 




*-T + i/ + -)' (96) 

where the ± symbol is changed to the — sign so that 6N remains a perturbation, and the 
trajectory o = is chosen. Coming back to the variables 5<p* and 5cr+ we see that Eq. (|96l) 
reproduces the x-dependent part of Eq. ([291) in view of Eqs. (1501) and (15T|) up to second order 
in 5d>+ and 5<r+: 




SN^^- _p + " exp[2JV(^-^)] + .... (97) 



Eq. (J96J), although reliable only up to second order, tells us that the expected behaviour of 5N 
in terms of 77^, 770-, x, and 7/ is indeed obtained. Moreover, from our previous discussion we know 
that 5N can be exactly written in terms of a series of x and y withouth the explicit appearance 
of the concerning exponential factors containing N. This is indeed partially confirmed up to 
third order when introducing Eqs. (1501) . (|5T|) . and (|52|) into the x-dependent part of Eq. (|29|) : 

^V + J) . (98) 




The bottom line of this discussion is that we have been able to identify two quantities that 
determine the truncation of the series up to some specific order. These two quantities are x and 
y which we could identify as the "coupling constants" of the theory in the context of Quantum 
Field Theory. By making \x\ <C 1 and \y\ <C 1 we can see from Eq. ( 1981) that all the terms higher 
than second order in x and y are subleading compared to the second-order ones. As regards the 
first-order terms compared to the second-order ones, we see that the latter are not necessarily 
subleading compared to the former because of the non-existence of the first-order y term and 
in view of \y/x\ < 1600 from Eqs. ( 1891) and (1901) and the values for 77^, 7] a and N considered in 
Sections [5] and El In the language of the Feynman-like diagrams [54J, truncating 5N in Eq. (1981) 
up to second order in x and y means considering only the leading diagrams at tree level and one 
loop which is what we have done in Sections [5] and [61 In fact, \x\ -C 1 and \y\ <C 1 mean that 

\x\ . ^«(£U<1, (99) 




1 1/2 f r? 3 (H \ 2 1 1 1/2 
■ exp[2iV(^ - Va )} « if M - exp[2iV(^ - Va )} « 1 , (100) 



2tt 



which are well satisfied for the cases when is dominated by the tree-level term (see Subsub- 
sectionE221- Eq. ([69]) and Subsubsection [5X2] - Eq. ([M])): 



;§)i-H^ i/2 -io-, (101) 

|(§) 2 ^ eXp[2iV( ^~^ )] } «{| ex P[-2^(l^|-|^l)]} <2- (102) 
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By explicitly calculating the two-loop and three-loop diagrams for and B^ and employing 
the results of Ref. [88], we have checked that the conditions \x\ <C 1 and \y\ -C 1 efectively make 
these diagrams subleading compared to the leading ones at one-loop level. 

Finally, we will discuss the convergence of the ( series in view of Eqs. ( 1931) . (Ell), and ( 195*1) . 
We first note that the series in Eq. ( 1941) is always convergent. As regards the series in Eq. ( 1951) . 
it will not be convergent at all while the function 

(103) 



{l + x) 2 + y 2 {l + x) 2 + B 2 x 2 ' 
with 

B= ^ expfiV^-^)], (104) 



does not satisfy the following necessary condition [87]: for the Taylor series around x = of a 
function f(x) to be convergent, it is necessary that the extension f(z) to the complex plane of 
f(x) is continous in a neighbourhood of z = 0. If this is the case, and the Taylor series of f(z) 
is indeed convergent, the convergence circle must either match or be inside the aforementioned 
neighbourhood. Of course, this is not a sufficient condition, but at least gives us a constraint 
on the possible values that x may take. 

Applying this condition to the expression in Eq. (11031) . we see that the extension of this 
function to the complex plane has poles for (1 + z) 2 = —B 2 z 2 which leads to 

±iB - 1 . . 

* = -p+r- (105) 

Therefore, the extension to the complex plane of Eq. (I103[) is continous for 



so the necessary condition for the convergence of the series in Eq. (195]) . and therefore for the 
convergence of the series in Eq. ( |93l) which is what we are interested in, is given by \x\ < 1. 
Thus, such a necessary condition for the convergence of the ( series is automatically satisfied 
once we choose |x| <C 1, as we have seen above it is required for working in a perturbative 
regime. 



8 Conclusions 

Observational cosmology is in its golden age: current satellite and balloon experiments are 
working extremely well 0H], dramatically improving the quality of data [6j. Moreover, foreseen 
experiments [7J [8] will take the field to a state of unprecedent precission where theoretical 
models will be subjected to the most demanding tests. Given such a state of affairs, it is 
essential to study the higher order statistical descriptors for cosmological quantities such as the 
primordial curvature perturbation (, which give us information about the non-gaussianity in 
their corresponding probability distribution functions. 
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( and its associated non-gaussianity depend on the specific inflationary model that describes 
the dynamics of the early Universe, the slow-roll class of inflationary models with canonical 
kinetic terms being the most popular and studied to date. Inflationary models of the slow-roll 
variety predict very well the spectral index in the spectrum P^ of ( but, if the kinetic terms are 
canonical, they seem to generate unobservable levels of non-gaussianity in the bispectrum 
and the trispectrum of C, making them impossible to test against the astonishing forthcoming 
data. Where does this conclusion come from? The answer relies on careful calculations of the 
levels of non-gaussianity f^L an d t nl by making use of the SN formalism [HJ HI [201 [22]. In 
this framework, ( is given in terms of the perturbation 5N in the amount of expansion from the 
time the cosmologically relevant scales exit the horizon until the time at which one wishes to 
calculate (. 

Due to the functional dependence of the amount of expansion, ( is usually Taylor-expanded 
(see Eq. ( 1291 ) and truncated up to some desired order so that f NL and r NL are easily calculated 
(see for instance Eq. (1551) ). Two key questions arise when noting that it is impossible to extract 
general and useful information from the ( series expansion in Eq. until one chooses a 

definite inflationary model and calculates explicitly the N derivatives. First of all, when writing 
a general expression for f^L or t nl i n terms of the N derivatives, how do we know that such 
an expression is correct if the series convergence has not been examined? Moreover, if the 
convergence radius of the £ series is already known, why is each term is the ( series supposed 
to be smaller than the previous one so that cutting the series at any desired order is thought to 
be enough to keep the leading terms? Nobody seems to have formulated these questions before 
and, by following a naive line of thinking, Jnl and t^l were calculated for slow-roll inflationary 
models with canonical kinetic terms without checking the ( series convergence and keeping only 
the presumably leading tree-level terms [T8] fl9l 120] 1221 [271 128] . 

These two questions have been addressed in this paper by paying attention to a particular 
quadratic small-field slow-roll model of inflation with two components and canonical kinetic 
terms (see Eq. (1431) ). Although the non-diagrammatic approach followed in Section [7] to find 
the necessary condition for the convergence of the ( series in our model might not be applicable 
to all the cases, we have been able to show that not being careful enough when choosing the 
right available parameter space could make the ( series, and therefore the calculation of /jvl 
and tnl from the truncated series (e.g. Eq. (13"3"j) ). meaningless. We also have been able to show 
in our model that the one-loop terms in the spectrum P^ and/or the bispectrum of ( could 
be bigger or lower than the corresponding tree-level terms, but are always much bigger than 
the corresponding terms whose order is higher than the one-loop order. If both P^ and B^ are 
dominated by the one-loop terms, a huge Jnl is generated which overwhelms the observational 
constraint, ruling out the model by an excess and not by a shortfall. If B^ is still dominated by 
the one-loop correction but P^ is now dominated by the tree-level term, sizeable and observable 
values for f^L are generated, so they can be tested against present and forthcoming observational 
data. Finally, if both P^ and B^ are dominated by the tree-level terms, J'nl is slow-roll suppressed 
as was originally predicted in Refs. [T51 EE ED] • 

What these results teach us is that the issue of the ( series convergence and loop corrections is 
essential for making correct predictions about the statistical descriptors of ( in the framework of 
the 5N formalism, and promising for finding high levels of non-gaussianity that can be compared 
with observations. In fact, now that we have learned the lesson, the level of non-gaussianity t^l 
for the same slow-roll model studied here will be the subject of a companion paper [35] . 
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Figure Al: Tree-level Feynman-like diagram for P^. The internal dashed line corresponds to a 
two-point correlator of field perturbations. 
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Appendix 

A Tree- level and one- loop diagrams for and 

We show in this appendix the mathematical expressions for the tree-level and one-loop Feynman- 
like diagrams associated with the spectrum P^ and the bispectrum of (, following the set of 
rules presented in Ref. [54]. To this end we have taken into account the N derivatives for our 
small-field slow- roll model given in Eqs. fl50l) . flBT]) . and (1521 . After presenting the mathematical 
expressions, we will estimate the order of magnitude of each diagram in order to determine the 
respective leading terms at tree-level and one- loop for both P^ and B^. 

Al Tree-level diagram for 

Looking at Fig. IA11 we see that Pf. ree is given by 

Pf ee = N*P„{k) 

2vr 2 1 /H^ 2 



fc 3 TlM V2tJ 



(Al) 

Of course, there is only one tree-level diagram for P^ and therefore Eq. flAlj) is the associated 
leading tree-level term. 

Our calculation in this appendix goes up to the one-loop diagrams so, in order to have 
complete consistency in the calculation [89], we should also take into account the one- loop 
correction to the two-point correlator in the field perturbations when calculating the diagram in 
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(a) 



(b) 



Figure A2: One-loop Feynman-like diagrams for P^. (a). The two internal dashed lines corre- 
spond to two-point correlators of field perturbations, (b). The internal dashed lines correspond 
to a three-point correlator of field perturbations. 



Fig. IA1I Such a correction has been studied in Refs. [90] ED [92J, [93l [91] where the most general 
result for single-field slow-roll inflation with N tota i not very much bigger than 62 is [91] 



P 



l—loop 



2vr2 



2n) \ + \2nm P 



-7T HkL) + 
o 



(A2) 



where L is the infrared cutoff for a minimal box [521 [53], and (5 is a renormalisation scheme- 
dependent constant that is expected to be negligible on large scales compared to ln(A;L) ~ 0{1). 
The one-loop correction to the field perturbation spectrum in Eq. (IA2[) is, therefore, negligible 
compared to the tree-level contribution P^ ee = (2n 2 /k 3 )(H+/27i) 2 if H+ <C m P as usually 
required. In our model if* <mp is indeed given but, since we are dealing with a two-component 
model, the result in Eq. (1A2j) may not be applicable. Anyway, we feel quite confident that the 
(up to now unknown) extension of Eq. (1A2j) to the multiple-field case will yield similar results, 
so we will keep the expression in Eq. (lAip as the leading tree-level contribution to P^. 



A2 One-loop diagrams for Pr 



Looking at Figs. IA2a and 

pi— loop b 



we see that P^ loor ; 



is given by two contributions P, 



l—loop a j 

and 



pi — loop a 



■pi— loop b 



1 r 
2 
1 

2 



NL + N; 



d 3 q 



P^(g)P^(|k+q|) 



+ 



1 

^1 



LVJ <i<>l 
d 3 q 



exp [4A^(?7 -r} a )] 



An 



d 3 q 



m X/h (k, q, |k + q|) + 
B s< j, 8a Sa(k,q, |k + q|) 



(27r) 3 ^ 5<t> 5<t> 
d 3 q 



(27T) 



4vr 



perm 



,1/2 



A4 3 (M,|k + q|) 



2tt/ 2 v / 2m P fc 3 g 3 |k + q| 



+ 



: U3 



exp [2^(77^ - r/ CT ) 



(A3) 



^ (2^) 3 pe ^ 2a .V2W 2V2m P fcV|k + q| 3 



M 3 (k,q, |k + q|) 



(A4) 
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where the ln(fcL) ~ factor comes from the evaluation of the momentum integrals in a 

minimal box [521 [53], the M.z{k\, k 2 , k 3 ) function is defined by 



M 3 (h, k 2 , h) = -hkl - 4^p + l -k\ + k 3 ) , (A5) 



with k t = ki + k 2 + k 3 , and the sub index perm. 12a. means a permutation over the last two 
arguments in .M3. 

A quick glance reveals that the first term in Eq. (IA3I) is subleading with respect to the 
second one because \r) a \ > |r/^| and exp[4Af(r^ — 77^)] ^> 1. The same is true for Eq. (IA4I) where 
exp [2N(rj^ — rj a )] ^> 1. Now, by comparing the orders of magnitude of the leading terms in Eqs. 
( 1A3I) and ( 1A4I) . we conclude that: 



pl-loo P a exp [ 4iV (^ _ Va) ] (B*y 



_c _ 

^ ^ex P [2iV(^-^)](f) ^ 

= ^^exp[2iV(^-r /(7 )]^»l, (A6) 
where mp> 0* and e^, < 1. Thus, the one-loop leading term for in our model is given by 

P V°° P = ^jLexp[4iV( W -^)] (AT) 



Having presented the leading tree-level and one-loop contributions to P^ in Eqs. (lAlj) and 
7]), a consistency issue to think about is the dependence of the expression in Eq. flA2j) on 



the infrared cutoff L. This quantity is in principle an artefact of the series expansion, and the 
final series result should in principle be independent on the chosen value for L (see for instance 
Ref. [95]). In fact, by assuming that this is the case, Refs. [52l 1551 [96] have shown that there 
is a running on the derivatives with respect to L so that changes in the ln(fcL) factors are 
compensated by the running of the A^ derivatives. This is similar to what happens in Quantum 
Field Theory where physical results independent on the energy scale must be independent of the 
chosen value for the renormalisation scale Q. Changing Q only modifies the relative weight of 
the tree-level and loop contributions, usually making the tree-level terms dominate over the loop 
corrections if Q is chosen around the relevant energy scale of the process studied. Nevertheless, 
we see that the ln(fcL) term in Eq. flA2j) does not compensate for the ln(fcL) term in Eq. 
(I A7I) . which is a real concern as we could expect since ( and its spectral functions are a set of 
observables. The solution to this paradox relies on the fact that the observed ( depends on L as 
the stochastic properties of the distributions depend on the size of the available region in which 
we are actually able to perform observations. In this regard ( is analogous to for instance the 
fine structure constant in Quantum Field Theory which, being an observable, depends on the 
energy scale for which experiments are done and, therefore, on Q. Likewise, ( and its spectral 
functions, though being observables, depend on the size of the regions where observations are 
done and, therefore, on L. Having this in mind it is essential to work in a minimal box [55J, i.e. 
with L a bit bigger than Hq 1 (with the subscript meaning today), so that ln(fcL) ~ 0(1) as 
has been done throughout this paper. 



29 



(a) (b) 

Figure A3: Tree-level Feynman-like diagrams for B^. (a). The two internal dashed lines corre- 
spond to two-point correlators of field perturbations, (b). The internal dashed lines correspond 
to a three-point correlator of field perturbations. 

A3 Tree-level diagrams for Bq 

Looking at Figs. IA3h and lA3b . we see that B l £ ee is given by two contributions B^ ree a and B^ ree b : 
B tree a = [p^{k x ) p^(fc 2 ) + 2 permutations] 

- -&ffifM§) 4 - (A8) 

B treeb = S^h, k 2 , k 3 ) 

Ul ±L\2nJ 2v / 2m P Uikf 



Now, from comparing the order of magnitude of the expressions in Eqs. (1A8|) and f]A9[) . we 
conclude that: 



B f eeb i Mv r^ 4 fT (Es 

= ( ai °) 

which in fact is usual as demonstrated in Refs. [T8J [97] . Thus, the tree-level leading term for 
B^ in our model is given by: 

M) • (Au) 

As was done for P^ in Subsection IA11 the one-loop correction to the spectrum of the field 
perturbations must be taken into account for the sake of consistency when calculating the 
contribution associated to the diagram in Fig. IA3h . The discussion about the relevance of 
this quantum one-loop correction is actually the same as in Subsection IA1I and, therefore, we 
may conclude with some confidence that the expression in Eq. (1A8I) is reliable. As regards 
the diagram in Fig. IA3b . it is necessary to include the one-loop correction the three-point 
correlator of the field perturbations in Eq. (1A9[) . which in fact nobody has calculated yet even 
for the single-field case. Nevertheless we might conjecture that, analogously to that for the P^ 
case, such a correction is negligible compared to the tree-level contribution to B^ and, therefore, 
the expression in Eq. (IA9I) will also be reliable. 
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(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



Figure A4: One-loop Feynman-like diagrams for B^. (a) and (b). The three internal dashed 
lines correspond to two-point correlators of field perturbations, (c), (d), and (e). The internal 
dashed lines correspond to a two-point and a three-point correlator of field perturbations, (f). 
The internal dashed lines correspond to a four-point correlator of field perturbations. 



A4 One-loop diagrams for 

Looking at Figs. IA4a . IA4b . IA4b . IA4d . IA4e . and IA4f . we see that B^ ree is given by six contribu- 
tions B\- loap \ Bj- ,opp6 , B\- loop \ B l - loop \ B\~ loap \ and B\- loap f : 



Tjl—loop a 



B 



l— loop b 

c 



pi— loop c 



d 3 q 
(2^ 



■P^^P^dki + qDP^dka-ql) 



' +4Lexp[6iV(^-^)] 



3^6 



gJ) MU0 (§)V (A12) 



- [N^N H N H4> + N^N^N^] x 



x 

1 

2 



d 3 q 



P Sf p{q)Ps^(\k 3 - q\)P s< j > {k 2 ) + 5 permutations 



(2tt) 3 
2 2?7 3 

d 3 q 



(A13) 



(2vr) : 



(q, |k 3 + q|, /c 3 )P < 5 9i (|ki - q|) + 2 permutations] + 
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J (2tt) 3 ^ p ^ m \2n) 2 v / 2m P g 3 |k 3 + q| 3 A;f |k x - q| 3 
+2 permutations J + 



4 ^ 



x 



7 — exp [47V(^ - ?7 CT )] x 

V / perm. 12a. 



e/ 2 M 3 (k 3 ,q, 


k 3 4 


-q|) 


1 


2 v / 2m P A;|g 3 


k 3 


+ q 


3 


ki - q| 3 



+ 



+2 permutationsj . (A 14) 

1 /" c? 3 ^ 

2 N 4> N H J Jc^y ^ Bs<p H s ^ k? " q ' ' ks ~ q D jP ^( /i;2 ) + 5 P erm utations] + 

1 f cftq 

+ -N lf) N ( f )( f > N aa J [Bs<t> 8a Sa{h, q, |k 3 - qD-P,^^) + 5 permutations] 

1 r f d*q 6 v (H*\ 6 e l J 2 M 3 (k 3 , q, |k 3 - q|) 1 
2rjffi U (2vr) 3 n p ^ m \2n) 2V2m P fc 3 V|k 3 - q| 3 fcf 

+5 permutationsj — 



/#*\ 6 A^3(fc3,g,|k 3 -q|) 1 | 

™ ,o„ V2tt/ 2A/2mp A;?o 3 |k 3 -q| 3 fc 3 



(2^) 3 * ^a^J 272m P fc 3 V|k 3 -q| 3 A; 3 
+5 permutationsj . (A15) 

-N^N^ J 3 5 ^ ^(fci, g, |ki + q|)P 5 4,(fc 2 ) + 5 permutations] + 
1 z 1 c? 3 g 

+-N%N aa(j) J -j—^ [B s< j> 8a 8a(h, q, |ki + <\\)P &(j> (k 2 ) + 5 permutations] 

1 r f d 3 q 6 /ff^ 6 6^ 2 -MaCfci.gJ^ + ql) 1 
(2vr) 3 n p ^ m \2nJ 2V2m P A; 3 g 3 |k 1 + q| 3 k\ 

+5 permutationsj — 
-^exp[2JV(ifo-ifc)]x 

_^_ 8 6 V f^V e ^ • M 3(fci,g,|k 1 + g|) 1 

(2vr) 3 * pJha^J 2v/2mp A; 3 g 3 |k 1 + q| 3 k\ 



U3 + 

ri-o 



+5 permutationsj . (A16) 
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d 3 q 



-I — Mi N 
2 ./ (2tt) 3 
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1 



■57T 
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, 4 exp [2N(ij^ - r/ ff )] x 
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+2 permutations 



277. 
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d 3 q 



6 V /^A 6 ^4(ki,k 2 ,q,k 3 -q) 1 

i^r 71 *<L 12a. V2 J A# 2 3 g 3 |k 3 - q| 3 ml + 



perm. /2a. 12a 

+2 permutations 



(A17) 

where the subindex perm. f2a. Via. means a permutation over the first two arguments and 
simultaneously over the last two arguments in M.±{k.\, k 2 , k 3 , k 4 ) defined by 



_M 4 (ki,k 2 ,k 3 ,k 4 ) 



^12^34 kt 



Ikl 



Zx2 ■ Z34 3 
— h 2k 2 • Z 34 + -(T12CT34 



2H 



034 



31 



"34 

W 124 4 



ki • k 2TIA , A;f 



with ky = kj + kj, k t = ki + k 2 + k 3 + k 4 , and 

,2 



a 



kj • kj + fcj 



— \~ 2 A?? /c a 

1 + — - + 3 



kl ' 



W, 



Iran 



j + & m "I - _|_ l(kik m + /c;/c n + k m k n ) ^ 6k[k m k n 



(A18) 

(A19) 
(A20) 

(A21) 
(A22) 



fe? fef 

Following the same kind of analysis as we carried out for the one-loop diagrams of and 
the tree-level terms for we conclude the following: 



j^l—loop a 
jjl—loop b 



Ei>€\ (h 



^ex P [4iV(^-^)]4^ (gjf) (§)' 
exp [2JV(^ - tfc)] > 1 , 



(A23) 



Tjl—loop a 



r>l — loop c 



j|,exp[6iVfa-^)]^)(f) 4ir* 
^ exp [4tffo - 4x* (f ) 6 g 
— ^ exp [2iV(^ - 7^)] > 1 , 



(A24) 
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rjl-loopd ~ \6e 1/2 /V fr 3 



^ ^£exp[4iV(^-^)]i»l, (A25) 



rjl-loope ~ o /rT x6 e 1/2 /V fr 3 



e 



— — exp [4JV(7fo - 77,)] > 1 , (A26) 



^ jlexptGiVfa-^)] (gig) (f) 6 4^ 



' \ 2 / \ 2 



exp [4A^ - r^)] > 1 . (A27) 



\V<pJ \ <P* 

Thus, the one-loop leading term for B^ in our model is given by: 

BV°" = jgj «P - *)1 (§)' K")^ (ff) ■ (A28) 



Once again, the ln(fcL) dependence in Eq. (1A28|) does not look like that obtained from 



introducing Eq. (1A2I) into Eq. (1A11I) . However the situation here is the same as that discussed 
at the end of Subsection IA2| leading us to identical conclusions. 
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